A statistical physics of stationary and metastable states: 
description of the plasma column experimental data 



A. Caboi, S. CurileP, A. Gonzalezi, N. G. Cabo-Bizet^.s and C. A. Vera^ 
^ Theoretical Physics Department, Institute de Cibemetica, 
Matemdtica y Fisica, Calle E, No. 309, Vedado, La Habana, Cuba. 
^ Departamento de Fisica, Universidad Catolica del Norte, Av. Angamos 0610, Antofagasta, Chile. 
- - - Departamento de Fisica, Centro de Aplicaciones Tecnologicas y Desarrollo Nuclear (CEADEN), 

Calle 30, esg. a 5ta Ave, Miramar, La Habana, Cuba. 
"^Instituto de Fisica, Universidad de Antioquia, Medellin, Colombia. 
^ Physikaliches Institut der Universitdt Bonn. Nussallee 12, Bonn, Deutschland. 



o 
o 



We propose a statistical mechanics for a general class of stationary and metastable equilibrium 
^ , states. For this purpose, the Gibbs extremal conditions are slightly modified in order to be applied 

• to a wide class of non-equilibrium states. As usual, it is assumed that the system maximizes the 

' entropy functional S, subjected to the standard conditions; i.e., constant energy and normalization 

CO , of the probability distribution. However, an extra conserved constraint function F is also assumed 

to exist, which forces the system to remain in the metastable configuration. Further, after assuming 
additivity for two quasi-independent subsystems, and that the new constraint commutes with density 
matrix p, it is argued that F should be an homogeneous function of the density matrix, at least 
^ . for systems in which the spectrum is sufficiently dense to be considered as continuous. The explicit 

^ ' form of F turns to be F{pi) = p', where pi are the eigenvalues of the density matrix and g is a 

a I real number to be determined. This q number appears as a kind of Tsallis parameter having the 

interpretation of the order of homogeneity of the constraint F. The procedure is applied to describe 
the results of the plasma experiment of Huang and DriscoU. The experimentally measured density 
is predicted with a similar precision as it is done with the use of the extremum of the enstrophy 
C/3 , and Tsallis procedures. However, the present results define the density at all the radial positions. 

In particular, the smooth tail shown by the experimental distribution turns to be predicted by the 
procedure. In this way, the scheme avoids the non-analyticity of the density profile at large distances 
arising in both of the mentioned alternative procedures. 
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I. INTRODUCTION 



^ ^ Nowadays, a relevant topic of research in statistical physics is the characterization of metastable and non-equilibrium 
^ . states of physical systems [l|, 0, H, 0, [1]. Multiple efforts have been made in order to explain this kind of behavior, 
' but the discussion intensively continues because of the large variety and complexity of physical processes involved 

m : |,i0,i. 

OO ' In this work we intend to present a more complete discussion of the statistical description previously sketched in 
Ref. [9]. We expect to describe the statistical physics of a large class of metastable and stationary states. Our 
00 ' main suggestion to deal with these situations is that statistical properties of the metastable or stationary non- 
' equilibrium states depart from the Gibbs thermal states in equilibrium in the following sense: the metastable or 
non- equilibrium configuration may be the results of an additional dynamical constraint, maintaining the system for 
a while in the metastable state, if the restriction is only an approximate one; and retaining the system in a time 
independent configuration, in the case of non- equilibrium stationary states. After the removal of the condition, the 
system is allowed to evolve in a thermal state described by the Bloch density matrix. Hence, the time independence 
of the statistical properties of the systems leads to commutation between the density matrix and the Hamiltonian 
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during the time for which the metastable or stationary state exists. As it was mentioned, we will assume that an 

extra constraint F is conserved in time, then F will also commute with the Hamiltonian. In the large class of systems 
having a non-degenerate spectrum of the energy, the Hamiltonian H , the density matrix p and the extra constraint 
F can be diagonalized in the same common basis of eigenstates. The same property is not necessarily valid when the 
energy spectrum is non degenerated. From now on, we will simply restrict the discussion to systems exhibiting the 
mentioned commutativity among H , p and F . 

Let us assume that the extra constraint defines the metastable state. Then, it is possible to note the existence of 
some related conserved constraints. A particular interesting one is given by the product operator HF . Note, that 
the assumed commutativity of p with F implies the possibility of expressing as a function of p. The above remark, 
motivated in us the idea of constructing a special analytical form of the condition to be incorporated in the Lagrange 
multiplier scheme for the maximization of the entropy. This condition will be constructed as a modified expectation 
value of the energy, in which in the density matrix in the usual mean value is substituted by the constraint function 
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F. This quantity results to be conserved as a direct consequence of the vahdity of the constraint F. 

As a result of the above construction, it follows that, when the additivity condition of the resulting statistical 
description is assumed to be valid for two approximately independent subsystems, the constraint function F should 
have the Tsallis homogeneous structure F{p) = Cqp'^ where g is a real number to be determined. It should be 
underlined that this conclusion follows for systems in which the energy spectrum is sufficiently dense for to be 
considered as continuous. For such systems, the analysis, seems to indicate an interpretation of the Tsallis q parameter 
as corresponding to the degree of homogeneity of the constraint F after to be represented as a function of the density 
matrix p. 

In order to start investigating its implications, the procedure was applied in this work to describe the experimental 
data obtained from the plasma experiment performed by Huang and DriscoU Q . The measurements of the electron 
density were chosen as an input of an iterative procedure for solving the equations for the probability density p. The 
results predicted for the densities at radial distances at which the their values are not small, describe the experimental 
results with similar quality as the ones following from two theoretical schemes existing for this problem. They are 
the minimization of the enstrophy on one case and the maximization of the Tsallis entropy in the other one p, 0] . 
It is known that these procedures are equivalent and their applications to the considered problem are recognized 
as main studies existing in recent literature 7]. In connection with the zone of small radial distances, both of the 
analysis describe this region with similar precision and their results are close to the experimental data. However, 
at large radial distances, where the density is small enough, their results strongly deviate from the experimental 
measurements reported by Huang and DriscoU [6(] . Our results radically improve the outcome of the aforementioned 
theoretical procedures in this large distance region. In this sense, the presented approach gives a new description of 
the smooth vanishing tail of the density profile, which has been experimentally measured. Each one of the alternative 
procedures needs to justify the presence of a non-analytical behavior of the density distributions in the region of the 
tail of small densities, as coming form the requirement of the non-negative character of the electron density. In the 
here proposed discussion, such a procedure is not becoming necessary and, as noted before, the dependence of the 
density as a function of the radial distance is predicted to smoothly decay. 

Our procedure has been applied to the quasiequilibrium state of excitonic polaritons ,10]. These polaritons are 
quasibosonic quasiparticles arisin g fr om the strong interaction of excitons and confined light modes in a microcavity. 
A finite model is studied in Ref. [lO| with a help of a master equation for the density matrix. The result is that the 
density matrix approximately commutes wit the Hamiltonian. We assumed that the quasiequilibrium is the result 
of a constraint in phase space and showed that an analytical expression for the density matrix, resulting from our 
procedure gives a reasonable fit to the numerical results. 

The paper proceeds as follows. In Section II we present the basic elements of the proposed statistical description. 
Section IIIIl then continues by arguing that the assumption of additivity in the statistical description implies that 
the constraint function should be a simple q power of the density matrix. In Section IIVI the procedure is applied to 
describe the results of the plasma experiment of Huang and DriscoU. Finally, Section |V] resume the results of the 
work. 



II. A STATISTICAL MECHANICS FOR A CLASS OF METASTABLE AND STATIONARY STATES 



Let us consider a physical system having a quantum dynamics described by a Hamiltonian H. In the Gibbs 
approach, the properties of the system in thermal equilibrium are contained in the Bloch density matrix 

which satisfies [i/, p] = and determines the conditional maximum of the entropy functional 

S = -Tr[p\og{p)] + a{Tr[pH] - E)+ 
+ /3(Tr[p]-l), 

where a and (3 are Lagrange multipliers corresponding to imposing two conditions: the conservation of the energy 
E and the normalization of the trace of the density operator p. The latter condition is placed in order to furnish a 
probability interpretation to the diagonal elements of p, this is Tr[p\ = 1. Eventually, additional conserved quantities 
can be added by also introducing their corresponding multipliers. 

Our main assumption in this work is that during a large relaxation time r, the considered physical systems are 
forbidden to reach the thermal equilibrium state, as a consequence of the action of extra constraints. These restrictions 
are assumed to be dynamically generated and able to obstructs during macroscopic time intervals the standard 
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evolution which normally drops the density matrix to the Bloch form. In other words, we propose the existence of a 
conserved quantity F, whose effects over the motion of the system is to delay the usual evolution, then leading the 
system to an intermediate equilibrium state differing from the Gibbs thermal one: the metastable or non equilibrium 
stationary state. Therefore, the system is compelled to remain in those states during the time lapse in that the extra 
constraint works. Thus, as a result of the precedent assumptions we will have: 

[H,F\ = Q. 

Therefore, the constraint F and the density matrix p both commute with the Hamiltonian. As mentioned before, in 
the large class of systems in which H has a non degenerated spectrum, it is possible to simultaneously diagonalizc H, 
p and F, within a common basis of eigcnfunctions. However, in a more general situation, we will simply restrict the 
discussion to the cases in that all the three quantities commute among them. Thus, the conserved constraint F can 
be expressed as a certain function of the density matrix F = F{p). Prom the mentioned assumptions if follows that 
the quantity F{p)H is also conserved 

[F{p)H,H]=Q. 

Further, the evolution of the system in the considered metastable or stationary states, also implies the time invariance 
of the specially constructed expectation value 

Tr[F{p)H] 
Tr[F{p)] ■ 

Now, lot us precisely state our main dynamical principle: the metastable and stationary states are detcrmincid by the 
maximization of entropy S subjected to the conditions of conservation of normalization, energy and the aforementioned 
expectation value. Then, the usual constraints plus the extra one, after multiplied by their corresponding Lagrange 
multipliers, are added to the entropy functional to construct the modified form of the Gibbs Lagrange multiplier 
scheme. The functional takes the explicit form 

S = -Tr[p log(p)] + a{Tr[pH] - E)+ 

= -Y,P^ logfe) + aiY,P^i^ ' E) + /3(^P^ - 1) + 
i i i 

where Pi,i = 1,2,3..., are the diagonal elements of the density matrix p = >< i\. 



III. TSALLIS q PARAMETER FROM THE ADDITIVITY OF THE DESCRIPTION 

Let us now consider implications of the modified statistical physics being introduced. They follow after assuming 
the additivity in the statistical description for the combination of two quasi independent systems. That is, each one 
of the subsystems will be assumed to be in the same kind of metastable state, showing analogous statistical properties 
as the whole body. Then, let us consider a pair of such systems that weakly interact between them. Any one of the 
bodies is considered in the same kind of metastable state. The entropy and the constraint functions of the first system 
can be written as 

-^pfhog(pf)) 
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Analogously, the same quantities for the second system may be expressed as 

i i 

,(2) _ L, ^U'. )f 



E 



Now let us consider that the description is also valid for the combination of both systems. Then, the entropy and 
constraints can be expressed as follows 

£(1,2) _ J1.2) (1,2) ^ _ (1,2) 



E 



^ (l,2)^ (1,2) 

(1,2) ^ 2^(1.]) ^ yP(i,l) 

i:(^.)^(p(S) ^ 



where the pair indicates the state of the composite system in terms of the indices of the states of the single 

systems as |i, j) — \i) x The assumed separability of the two bodies allows to write a relation for the probability 
of the combined states in the following way 

(1,2) (1) (2) 

= Pi P) 

T.^^ =i:p^t.p^ 

= 1. 

Thus, the probability constraint for the individual subsystems is valid for the probability constraint of the composite 
system. 

The additivity of entropies follows in a similar way 

ihj) 

i 

= 5(1) +5(2). 

We have the same outcome for the addition of energies 

B(1,2) _ (1,2) (1,2) 
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(l)J2)Jl)_^X^ Jl)J2)J2) 

(»,i) («,i) 

(2) (2) 



J 



= i;(i)+£;(2). 
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Let us come to the main new element considered in this section. Now we assume that the g-expectation value of the 
energy also satisfies the statistical independence condition for the combined system as 

F{p^^']{) = Fip'^^)F{pf). (1) 

Therefore, the g-expectation value also satisfies the additivity properties, as follows 



(1,2)^(1,2) 



F — ^ ^(1,2) ' 

[id) 

p(l)F(2)r^(i) I ^(2)^ 



Y- p(i)p(2) 
Y- p(i)p(2)^(i) Y" 



Y- p(i)p(2) Y" p(i)p(2)^(i) 

2^{t.3)^i ^3 ^(i,3)^^ ^3 



However, this condition also imposes a strong restriction on the possible forms of the function F defining the newly 
introduced expectation value. In order to see this central point of the presentation, the statistical independence 
condition ^ will be rewritten. Consider F as expanded in powers of its argument x, in the general form 



oo 



F{x)^x-Y.f^^''- (2) 



n=0 



Combining Eqs.|T]) and ([2]), it is possible to write 



oo oo oo 

71—0 Ti— m— 



After a suitable rearrangement of the variables, it follows 

= E E /"( Pfn^nm - frn). 

n— m— 

Next, let us assume that the domain of the probability pi for all the values indices i is a continuous set. Then, the 
completeness of the basis formed by the powers of a variable implies 

fniSnm " fm) = 0, for all m and n . 

At this point, let us suppose that fm 7^ 0, for a particular value of m = mo- Then, the validity of the above relation 
directly implies that fm = for all the values of m 7^ mo- Therefore, the explicit form of the operator as a function 
of p can be written in the following simple way 

= f, P' (3) 
= ^Pl\i><3\ 

i 

The parameter g is a real number which value should be dynamically determined. The above conclusion complete the 
argue indicating the Tsallis mean value structure of the modified average imposing the dynamical constraint F. 
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A. Extremal entropy equations 



Let us examine in this subsection the F dependent modified expectation value. As noted before, its most interesting 
property is the full coincidence with the Tsallis g-expectation value of the energy 

T7 — T7' ^iPi^i 

where the index q is the power of the probability in ([3]). The modified entropy functional takes the form 

S - -^paog(p,) + - E)+(3(^p, - 1)+ 



and the Lagrange extremum equations following from the functional are 

— ,* = 1,2,... = 0, 

opi oa 

d(3 ' 97 



Their explicit calculation leads to the following set of coupled equations for the eigenvalues of the density matrix p 



l-ae,-l3^~log{p,)+j^^ie,-Eq),i^l,2,... 



E = ^Pi<ii, 

i 

1 = ^K, 

i 

After some algebra, we obtain the integral relation 

= - Y,P^ logbO + -Eg)-l + aE + P, 

= ~^p, log(pi) -1 + aE + p 

i 

= S + aE -l + p. 

The last expression establishes a connection among the entropy, energy and the Lagrange multipliers. As usual, it 
seems useful to be employed in constructing generalizations of the free energy and other thermodynamical potentials 
and relations. However, we will not deal with this discussion here. 



IV. METAEQUILIBRIUM STATES IN ELECTRON PLASMA COLUMNS 

In this section, we apply the statistical procedure being investigated to an important example: the relaxation of 
a 2D turbulence to a metaequilibrium state. A relevant experiment related with this problem was performed and 
discussed by Huang and DriscoU [6]. They compared results from several theoretical approaches to their measured 
data. The theoretical scheme which furnish the best approximation, is the minimum enstrophy model. The experiment 
considered the 2D dynamics of a electron plasma fluid. The authors were able to identify and measure properties of a 
metaequilibrium state (MES) for a plasma sample in which a turbulence initially occurs within a rotating magnetized 
cylindrical electron column. The turbulence relaxes to a long-lasting metaequilibria state with axial symmetry. The 
relevant conserved quantities of the electron column in the MES are the following: energy H, entropy S, angular 
momentum Pg and number of particles N^. 
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1 R 

We start from the Poisson equation in order to relate the potential — to the particle density p = nj^, 

where cf) and n are the potential and the particle density (in the international system of units). The parameter e is 
the electron charge and the radius of the rotating drum. In agreement with Rcf. 6] all magnitudes are expressed 
in a special set of units, in which the cylinder radius is equal to 1. The Poisson equation for 5" is 

y2 ^ 47rp(r) (4) 
I d f d^(r) 



r dr \ dr 

The solution of this equation can be obtained in the form 

*(r) [ 2nr'dr'F^{r,r')p{r'), (5) 
Jo 

F(f(r,r') = 21n(r)e(r-r') + 21n(r')e(r'-r), 

where the gauge freedom of the potential has been employed to obtain a symmetric expression of the Green function 
i^if (r, r'). Note that this fact may introduce a change in the value of the total electrostatic energy with respect to the 
treatment in Ref. 6j. 

The entropy S, normalization N, energy H, angular momentum Pg and the constraint Eg are given by the expres- 
sions 



1 



fl" = -- / 27rrdr*(r)p(r), (6) 



1 

d^r{l-r^)p{r), 

N = 1= [ 2TTrdrp{r), (7) 

JO 

S = — 2nrdrp{r)\n[p{r)], 



Jp^{r)(-^) 



d\ 



J pi{r)d'^r 

We should then find the extremum of the functional 

p{r) ln[p{r)]d\ + as {H - 

+ aN (N - I) + ap {Pg - Pj^) 
+ l{E,- El) . 

The Lagrange multipliers are given by aE^ ctp, ajv and 7. The corresponding Euler equations are 

ln[pH] = {aN - 1) + ap{l - LU^) (8) 
- OiE j d^rF^{r,Lu)p{r) 



/ d'^rpi{r) V 2 
7 f 1 



-qp'^-'icj)^icj) 



d^rp'i{r)Fi{r,uj) 



J d'^rpi{r) V 2 

The integrated version of the above relation is 

S+{aN-l) + 2aEH + -fEg + apPg:=:0. (9) 
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We start from the measured data, for the electron density, given in Ref. and solve iteratively for p in a 
discrete version of relation ([5]) 

ln[p(wj)] = (aw-l) + «p(l-^') (10) 
- aE 27TriAriFQ {ri,ujj)p{ri) 
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^,27rrfcArfcp9(rfe) 

7 

J2t 27rrfcArfcp9(rfc) 

/ 1 NT^ 



2 ^ 27rr,Ar,p«(r,)i^o (^«. ^j) 



7 



X;27rrfeArfep9(rfc) 



to obtain consecutive approaches to the density distribution by exponentiating the logarithm in the left hand side 
(l.h.s.) of the relation (jTU]). Now, the indices i,j=l,2,... indicate the increasing values of all the radial positions 
(or Wi) at which the electron density measurements were reported in Ref. @. The initial density function p^{ujj) 
was constructed by assigning the density value measured at each point wj in Ref. @. The distance elements Ar^ 
are defined as the differences Ark = ''''^^ 2'''°^^ ' ^ ^ 2...N — 1. But Ari — r2 — ri and Ar^ = — t^-i- Let 
us denote the l.h.s. of Eq. pU]) evaluated in the density values obtained in Ref. 6] as L°, and the right hand side 
(r.h.s.) as i?^. The procedure starts by evaluating the vector — i?^ for all its components j, in the initial density 
data p^{LOj). Then, this vector becomes a function of the Lagrange multipliers and the value of q. Next, the values 
of those parameters a%,a^p,a%, 7^^ and that minimize the quadratic difference X]j(^j ~ ^^j)^ were determined. 
Finalizing this first step, the set of obtained multipliers were substituted in the expression for i?° together with the 
initial density p"(wj), defining in this way as the l.h.s. of pO|) . The new density data to be employed in the next 
iteration p^(wj) were determined from p^{ujj) — exp[Lj]/ 27rrfe Ar^ exp[L^]) , where the denominator shows that 
the normalization is imposed. The new step continues by writing a new vector Lj — i?j in which the density ^^(wj) is 
substituted in i?j and the parameters are set free again. Then, the expression X)j(^] ^ ^})^ again minimized to 
get the new optimal multipliers and q. This iterative procedure was carried out thirty times and the resulting degree 
of convergence is expressed by the value Ej (-^f " Rf f = 2-58 x 10"^ in comparison with the initial evaluation 
giving a result Y.j{L°j - ^j)^ = 0-33. Another indicator of the convergence is Eq. In this case the l.h.s goes from 
-6742 in the first step, down to 0.1 in the last one. Thus, the numerical solution of the system was taken as the 30*'' 
iterative results for the density values and the set of parameters a^, q;|?, aj^, 7"^" and 

In Fig[T]the data of Huang-DriscoU [6| are depicted in comparison with the ones predicted by the restricted minimum 
enstrophy model and our iterative solution. The values of the parameters obtained from the iterative procedure are 
given by 

af = (3.9±0.3)10^ (11) 

= (1.8±0.1)10^ 
af = -4.627 ±0.001, 
7^0 === (3.9±0.3)10^ 



^30 



0.999993 ±0.000001. 



Both for the parameters and the physical magnitudes, the errors are given by the difference between the value of 
the quantity in the last iteration 30*'* and in the previous 29*'*. The values obtained for the physical magnitudes are 

S = 0.0716522 ±0.000001, (12) 
H = 0.476445 ± 0.000003, 
Pg = 0.866057 ±0.0000001, 

which are reasonably close to the experimental values Q given the natural presence of experimental errors and our 
discrete approximation of the integrals. The experimental values are S^'^p = 0.087 and Pg^^ = 0.861. 

It can be noted that the results for the Lagrange multipliers associated to the Energy constraint, the Tsallis 
constraint like condition an the number of particles requirement are relatively high quantities. This curious outcome 
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FIG. 1: Density of particles p(r) vs. the radial position. The squares show the experimental result p^{rj) of Ref. Q. The 
continuous curve illustrates the density obtained from the minimal enstrophy model. The diamonds indicate the density 
evaluated from the iterative solution of the maximum entropy principle subject to the constraint Eq-. p'^^\rj) . 



results to be compatible with the resulting value of q which is very close to the unit. In this case, the only way in 
which the finite differences between the Gibbs and the present resuh for the density could arise is due to those large 
values for the multipliers. This is so because because the additional constraint for such a value of q very close to the 
unit is very close to the Gibbs energy constraint. In order to rule out the existence of a possible instability of the 
results for the multipliers, as a function of the step in the iterative process, in Fig. [2] we plotted the multipliers for 
the Energy and Tsallis constraints versus the iteration number. As it can be seen, the results continuously vary and 
decrease in moduli with the increase of the iterative step. 




FIG. 2: Lagrange multipliers as and 7 as functions of the iteration number. In spite of being relatively large numbers, their 
behaviors do not exhibit instabilities in their values when varying the step. The observed property is also fulfilled for the 
multiplier of the number of particles ajv. 

In order to search for the sensibility of the results on the way of evaluating the integrations, the iterative method 
described was carry out as well with a different way to approximate the integrals. They were now calculated by taking 
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a trapezoidal approach: g{r)dr « J27=i Ar^ 9(ri)+g(rz+i) ^ ^j^j^ _ Vk+i — r^. The obtained results are basically 
the same as the ones showed here for: the density solution of Eg. pH)) . the magnitudes of the lagrange multipliers and 
the convergence of l.h.s from Eq.® to zero. 

The Fig. [1] indicates that our solution reasonably agrees with the experimental data presented in Ref. 6] . At the 
points with bigger values of the density the results are close to the ones predicted by restricted enstrophy calculation. 
However, comparing the curves, one can appreciate that the minimum enstrophy model has an abrupt ending for a 
given radius tq, which is a limitation of this approach. In this case, the density is written as an analytical expression 
given by p{r) = a{Jo{f3r) — Jo(/3ro)), for r < ro or vanishes otherwise. In contrast, the present discussion correctly 
predicts the experimental smooth decay of the density to zero at large distances. Thus, in that region, the proposed 
procedure is able to avoid the singularity predicted by the minimal enstrophy analysis. It can be concluded that the 
performed numerical study shows that the statistical mechanics proposal investigated in this work, reasonably well 
describes the metaequilibrium state found in Ref. Q. 




FIG. 3; Comparison of the results for the density obtained by disregarding or not the Tsallis like constraint. The open squares 
show the Gibbs results, that is the solution of the iterative procedure by only retaining the fixed energy constraint. The 
diamonds again indicate the density evaluated from the iterative solution of the maximum entropy principle subject to both: 
fixed energy condition and the new constraint Eq. 

Finally, in Fig. [3]and for comparison purposes, the obtained density is plotted in common with the density evaluated 
from the same procedure but in which the Tsallis like constraint Eq is not considered: 7 = 0. That is: the conserved 
magnitudes are Energy, Angular Momentum and Normalization. Both results drastically deviate from each other 
at small radial distances. A recent work studying the behavior of a similar system [ll| establishes that in certain 
ranges of the ratio H/Pg it could dominate an evolution guided by a near maximum Entropy principle (for bigger 
ratio), or alternatively, a minimum Enstrophy (for smaller ratio). The Huang-DriscoU experiment furnishes a primer 
of the second case, since it is approximately well described by a minimum Enstrophy state. Therefore, as we are able 
to reasonable well describe their data our results deviate from the Gibbs behavior, and furnish support to the 
additional constraint argued here to be relevant in determining metastable and stationary states. 

It should be underlined that an issue that is yet puzzling for us, is the reason for that very small deviations from 
the unit of the q number (and its associated large values of the multipliers), are precisely the ones leading to a 
reasonable good match between the numerical solution and the quasiequilibrium plasma density. We expect to be 
able of considering this question elsewhere. 



V. SUMMARY 



In this work, we propose a modification of the Gibbs approach in order to deal with stationary and metastable 
equilibrium states of physical systems. As usual, the time evolution is assumed to maximize the value of the entropy 
S by satisfying the standard constraints, like constant energy and normalization. However, the description of the 
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no-equilibrium states is searched by considering that those configurations are stabihzed by the presence of an extra 
constraint F, being (approximately) conserved in the motion. The constraint, being restricted to be time invariant, 
is assumed to play the central role in stopping the evolution of the entropy to its maximal value in the Gibbs thermal 
equilibrium state. Afterwards, we assume a situation in which the description, which is valid for a composite system 
formed by two quasi independent subsystems, is also valid for each of them. This supposition implies that, if the 
constraint F commutes with density matrix, then F{pi) has a simplified dependence on pi of the form F{pi) = pi. 
One of the main outcomes of the work arises from this conclusion: The suggestion of an interpretation of the Tsallis 
q parameter as the order of the homogeneity of the constraint F, when it is expressed as a function of the density 
matrix. Remarkably, for small and medium radial distances the approach furnish results of similar quality as the ones 
given by the extremum of the enstrophy and Tsallis procedures. But, moreover, the smooth tail of the experimental 
density distribution at large distances is also predicted. Therefore, the analysis becomes able of avoiding the non- 
analyticity shown by the density curve in the mentioned alternative schemes, by describing the density decay tail of the 
Huang-DriscoU experimental data at large radial distances. Finally, it can be referred, that the proposed description 
has been also applied in Ref. [l3|, to a gas of polaritons contained within a quantum dot pumped by an optical 
resonator. Similarly as in the here discussed plasma problem, the scheme furnished a reasonably good description of 
the statistical distributions associated to that system. 
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